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We investigate the nonlinear dynamics of hybrid inflation models, which are characterized by two 
real scalar fields interacting quadratically. We start by solving numerically the coupled Klein-Gordon 
equations in static Minkowski spacetime, searching for possible coherent structures. We find long- 
lived, localized configurations, which we identify as a new kind of oscillon. We demonstrate that these 
two-field oscillons allow for "excited" states with much longer lifetimes than those found in previous 
studies of single-field oscillons. We then solve the coupled field equations in an expanding Friedmann- 
Robertson- Walker spacetime, finding that as the field responsible for inflating the Universe rolls down 
to oscillate about its minimum, it triggers the formation of long-lived two-field oscillons, which can 
contribute up to 20% of the total energy density of the Universe. We show that these oscillons 
emerge for a wide range of parameters consistent with WMAP 7-year data. These objects contain 
total energy of about 25 x 10^" GeV, localized in a region of approximate radius 6 x 10~^® cm. We 
argue that these structures could have played a key role during the reheating of the Universe. 

PACS numbers: 98.80.Cq, ll.lO.Lm, ll.15.Ex 



I. INTRODUCTION 



Inflation [T| is widely regarded as the mechanism behind the generation of the density fluctuations that eventually 
led to the formation of large scale structure in the Universe. In its simplest version, inflation is driven by a scalar 
field, the inflaton, which slowly rolls down its potential. While the total energy density is dominated by the inflaton's 
potential energy, the Universe grows exponentially fast. Eventually, the inflaton decays and reheats the Universe [2], 
giving rise to a radiation-dominated era of expansion consistent with big bang nucleosynthesis [3j . Given that we do 
not have a compelling candidate for the inflaton or for the fields it decays into during reheating, most of the work in 
inflationary cosmology focuses on the properties of different models and how they can be constrained by current data 

A particular class of slow-roll inflation models that is still viable by present constraints is provided by hybrid 
inflation [6j, where in addition to the inflaton 4>, there is another fleld ij), called the "waterfall" field. While the 
inflaton has a zero vev, the waterfall field has a nonzero vev. During inflation, the waterfall field is trapped at the 
minimum of its potential. After the inflaton crosses below a critical value (pc during its slow roll, the curvature at the 
origin becomes negative and the waterfall field undergoes spontaneous symmetry breaking: both fields are quickly 
driven towards their vevs, effectively terminating inflation and setting the initial conditions for reheating [7]. 

The dynamics of reheating after hybrid inflation has been studied both analytically and numerically [31 [SI E] . In 
this paper, we report results of 3d numerical simulations of the coupled scalar field equations just after the end of 
hybrid inflation, as the waterfall field ■0 relaxes spinodally towards its vev and the inflaton starts to oscillate about its 
minimum. Because of the quadratic coupling between the two fields, nonlinear effects are prevalent. Indeed, inflaton 
"hot spots" were reported in [9|, while "condensate lumps" were found for a class of supersymmetric models in [lOj . 
Here, we observe the emergence of long-lived, time-dependent field configurations, which we identify as a new class 
of two-field oscillons [Tmi4j . If oscillons are sufficiently long-lived, that is, if their lifetime is at least of order of , 
they can have important consequences for the cosmological dynamics. The potential role of oscillons in cosmology 
has been investigated in Id simulations [15H17j , while 3d simulations have been carried out in the case of a single real 
field in a double-well potential [18] and in the case of a 0® potential [19] where the emergence of flat-top oscillons 
[20] was studied. We find that these two-field, "hybrid oscillons" are indeed very long-lived, and emerge for a wide 
range of model parameters. Given that the inflaton must couple to other fields, and that oscillons have been found 
in a wide variety of models, from Abelian-Higgs models [21] to the standard model [22,, we should expect oscillonlike 
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structures to emerge during preheating, considerably enriching the post-inflationary dynamics and the reheating of 
the Universe. 

This paper is organized as follows. In the next section, we investigate the model with two coupled scalar fields, in 
the absence of cosmological expansion, and show that it admits spherically-symmetric oscillonlike solutions. These 
solutions are attractors in field configuration space: different initial conditions lead to long-lived, localized configura- 
tions with the same total energy. We also show that our model admits two different types of these oscillon solutions: 
the "ground state," where oscillons have similar lifetimes and energies to the well-known single- field models [T^HIl], 
and a newly found "excited state," where oscillon lifetimes can be three to four times as large as those for the ground 
state. We may think of the excited state as a metastable state with two possible decay channels: either it decays 
by directly radiating its energy to spatial infinity or, quite interestingly, it decays to the ground state configuration, 
where it remains for a while before radiating its energy to spatial infinity. This two-step decay inspires us to associate 
these higher-energy configurations with excited oscillon states. In Sec. |III|we introduce the hybrid inflation model and 



constrain the parameters to be consistent with WMAP 7- year data (23^. In Sec. IV we report results from numerical 
simulations in an expanding Universe, calculating the fraction of energy in oscillons for a range of values of the ratio 
of coupling constants g^/A, where is the quadratic coupling between the two scalar fields and A is the quartic 
coupling of the waterfall field. Moreover, we show that for g^/A = 2, oscillon production is particularly efficient, 
with the fraction of energy in oscillons reaching approximately 20%. We show that this effect is due to parametric 
resonance of a band of low k modes of the waterfall field, arising from the resonant oscillation of its zero mode and 
of the inflaton about its minimum. We conclude with a summary of our results and directions for future work. 

II. OSCILLONS IN A TWO SCALAR FIELD MODEL 

In [12H14j . spherically-symmetric oscillons were found in a single scalar field model with a symmetric (and also 
with asymmetric) double-well potential, starting from an initial Gaussian field profile that interpolated between the 
two minima of the potential. We perform a similar analysis in (3 H- 1) dimensions, now working with two scalar fields 
coupled quadratically, governed by a potential used in many hybrid infiation models 

= ^ - + ^^'-^^ + (1) 

Defining dimcnsionless variables 5;^ = Mx^, tp = ip{M / \/X)^^ , <}> — (j){M / , setting h — c — 1 and dropping the 
tildes, the equations of motion obeyed by the scalar fields become 

^-V^0 = -(^^0+^V^20) (2) 

72„;, _ fj, 9^ J,2„;, „/,3 



V, - = (^^ - y rV' - rj , (3) 

where an overdot denotes the derivative with respect to dimcnsionless time. For the moment, we will only consider 
spherically-symmetric configurations such that — ?► /dr^ + {2/r)d/dr in Eqs.[2|and [sj We look for oscillon solu- 
tions starting from Gaussian initial profiles for both fields, ^{r, 0) = 2cxp(— r^/i?^ — 1 and 0(r, 0) — 2 exp(— r^/i?^) 

while setting 0(r, 0) = '4){r, 0) — 0. The effective masses of the two fields (in units of M) are given by 

2 2 

We discretize the equations of motion using fourth-order spatial derivatives with lattice spacing Ar = 0.02 and 
step forward in time using a fourth-order Riingc-Kutta method with At = 0.01. We have verified that any further 
reduction of Ar and At does not affect our results. We add a damping term near the lattice boundary to absorb any 
radiation that may be reflected back and interfere with the oscillon evolution [53]. We ensure that no energy comes 
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back into the oscillon region by computing the flow of energy through the "energy shefl" and verifying that no energy 
goes back in. 

The equations of motion contain two free parameters, as of yet unconstrained: the ratio of the square of the two 
mass scales m^/Af^, and the ratio of couphngs Since the main goal of this paper is not to explore in detail the 

properties of two-field oscillons (a worthy task left for a future publication) but to show that they exist and explore 
their cosmological implications, in this section we fix the parameters as nn?/]\-P — 1.5 x 10^^ and g^/\ = 2.0, even 
though oscillons can be found for a wide range of these two parameters. The reason for this choice will become 
apparent in the next section. We note that with these values, the fundamental frequencies of oscillation near the 
minimum of the potential (at <?!> = and ^p = ±1) are almost equal, with = ^/2M and ~ ^/g^JXM = \/2M. 

We have solved the coupled nonlinear partial differential equations Eqs. [2]and[3]for a wide range of initial and 
the two free parameters in the initial Gaussian profiles of the fields, looKing for combinations that lead to oscillon 
formation. Once formed, an oscillon's energy is operationally defined as the energy contained inside a sphere of radius 
of i?sheU = 10M~^ around the origin 

There are four possible outcomes: 1. The fields do not settle into an oscillon state and quickly radiate their initial 
energy to spatial infinity; 2. The fields settle into an oscillon state and reach an energy plateau at E'osc ~ 43M/A, 
similarly to the single-field models studied in [T2II3]. These oscillons survive for lifetimes ranging from ~ 3, OOOAf 
to ^ 12,000M~^, depending on the initial Gaussian radii. We note that the longest-living single-field oscillon had a 
lifetime of ^ 8, 000M~^ [T31II1]. This oscillon state will henceforth be called the ground state; 3. The fields settle 
into an oscillon state with energy plateau E'osc ^ 69A//A and lifetimes up to ^ 33, OOOAf^^. This configuration will 
be called the excited state; 4. The fields settle into an excited oscillon state and then transition to the ground oscillon 
state where they remain for a while before decaying. In this case, we found that oscillons can remain localized for a 
total lifetime of ~ 40, 000M~^, some five times longer than their cousins in single-field models. The energy evolution 
of the three classes of oscillons is plotted in Figjl] 




FIG. 1: Oscillon energy versus time for three distinct classes of oscillons in our model. An oscillon in the ground state is plotted 
on the left. Two excited oscillon states are plotted on the right, one of which transitions to the ground state before decaying. 
Energy is given in units of A// A and time in units of . 

The two oscillon states can be contrasted by looking at the behavior of the fields at the origin {r — 0), as shown in 
Fig. [2j In the ground state, the two fields coherently oscillate about their respective minima, maintaining an almost 
Gaussian profile with ^ 3.25 and R^ ~ 3.2. The behavior of ip is similar to the one-field oscillon, leading to 
lifetimes enhanced by roughly ^ 30%. (Since we have not performed an extensive analysis for different values of the 
two parameters g^/X and m? /NP, these statements may be modified in the general case. It would be interesting to 
examine this question in more detail to see whether resonant behavior between the two fields may exist leading to 
much longer oscillon lifetimes.) 

In the excited oscillon state, the fields again reach a near-Gaussian profile with R^ ^ 2.0 and i?^ ~ 2.5. However, 
oscillates about the minimum of its effective potential with large amplitude, and so the g'^cjy^ijj'^ term in the potential of 
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FIG. 2: Oscillation of the fields i/)(r = 0, t), 0(r = 0,t) at the oscillon core, in the ground state (left) and in the excited state 
(right). 



Eq. [T| prevents ?A from reaching its minimum: the field V' remains "trapped" and periodically returns to its symmetric 
state at = 0, which is why we identify this configuration the excited (or metastable) state denomination. This 
behavior leads to considerably longer oscillon lifetimes, suggesting that if large-amplitude oscillations of the field ^ 
about its minimum could be maintained, oscillons may survive for arbitrarily long lifetimes. Figure [3] illustrates how 
the large-amplitude oscillations of affect the potential for ■0. The amplitude of required to restore the symmetry 
has to satisfy 



92y(<^,^) 



>0, 



(7) 



i/'=0 



which gives > \/A/.9 = a/2/2 for our choice of parameters. 




FIG. 3: Snapshots of the potential in the ■(/' direction for the ground state (left) and for the excited state (right) at different 
times during a period in an oscillon's life. In the excited oscillon case, the amplitude of is large enough (> \/2/2) to restore 
the symmetry and keep ^ "trapped" without reaching its broken symmetric vacuum expectation value. In contrast, during the 
lifetime of the ground state oscillon the symmetry is never fully restored. 

In Fig. |4) we show how the existence and lifetime of oscillons depends on the choice of the initial Gaussian radii. 
We vary both radii from 1.0 to 5.0 in steps of 0.05. The system portrays a subtly fine dependence on the initial 
conditions. Areas where the system does not settle into an oscillon are shown in black, whereas a nonblack color 
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indicates oscillons with different lifetimes. In Fig. [5] we separate the same data for ground state and excited state 
oscillons. 

For small radii, energy considerations alone clearly preclude oscillon formation. An obvious lower bound is 
E{R^,Rci,) > A3M/\, which lies in a band roughly given by + i?^ < (1.5)^. However, we note that no oscil- 
lons form for a band roughly bounded by R^ + R^ < (2.4)'^. A similar lower bound on the initial radius that gives rise 
to oscillons was obtained for a single oscillon case [13] ■ We can use the same approach to obtain an analogous result 
in the two-field model. First, write the fields as (j){r,t) = (/)o(i) exp[— r^/i?^] and ip{r,t) — ipo{t) exp[—r'^ / R'^] — 1. 
Integrating the Lagrangian density in the radial coordinate and varying the result with respect to the field ampli- 
tudes (j)o{t) and tpoit), we obtain two coupled ordinary nonlinear differential equations. Performing a small-amplitude 
expansion in both fields (</> = 0o + <50 a-nd "0 = "00 + ^ip) a-nd keeping terms linear in the fluctuations, we can show, 
focusing on Sip, that fluctuations will be unstable if 



3^2,2 4^6, 
— i-V'o 
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bl<o, 



(8) 



where 0o is the amplitude of the (j> field. These unstable fluctuations signal the existence of oscillons [TH[T3]. Given 
that is a parabola with minimum at % = 8v/3/9, we find that for oscillons to exist, we must have 
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< 0. 



(9) 



The inequality is thus consistent with the black band of Fig. [4j since the contribution from the field (p never acts to 
decrease the minimum oscillon radius at R^ = 2.42. For g"^ = we reproduce the result in Ref. [T^], R,p > 2.42. 
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FIG. 4: Oscillon lifetime for different combinations of initial and Rcj,. Radii within the black regions form no oscillons. A 
color gradient towards white indicates the range of growing oscillon lifetimes. 



III. HYBRID INFLATION: MODEL AND CONSTRAINTS 

We now investigate whether such two-field oscillons could play a role in cosmology. We turn our attention to the 
simplest realization of hybrid inflation |S] , governed by the potential described in Eq. [ij With a Friedmann- Robertson- 
Walker metric, the complete equations of motion for the two flelds become 



(10) 



..0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 5.0 



..0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 5.0 



FIG. 5: Breakdown of Fig. |4] with ground state oscillons on the left and excited state oscillons on the right. Although excited 
state oscillons are rarer, they live considerably longer. There are also clear linear relationships between the radii, suggesting 
resonance patterns that remain to be explored. 
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subject to the Friedmann equation. 
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(12) 



where rupi is the reduced Planck mass. 

For large values of 4>, the effective mass squared is positive for both (p and tp and the potential has the symmetry 
4> ^ —ip with a global minimum ai cj) = i/j = 0. For values of cj) below the critical point (j)c = M/g, the effective mass 
squared of tp becomes negative and symmetry breaking occurs, shifting the global minimum to ipc — M /\/\. Fast 
inflation is realized while (j) slowly rolls down its potential with ■0 = 0, and lasts until (j) drops below the critical point (j)c^ 
where symmetry breaking occurs and inflation ends. During inflation, the potential is T^((/), ip — Q) = M''/4A+m^(/)^/2. 
Since (j) is of the order 0c — M/g, for m? ^ g'^M^ /\ the vacuum energy dominates the total energy density. The 
expansion rate during inflation is then well approximated by Hf^^ ~ M'*/(12ATOp[) [71. In the limit m <C -ffinf, the 
amplitude of the curvature perturbation spectrum can be calculated exactly in the case of hybrid inflation j25j 



1/2 

n 



(13) 



and WMAP 7-year data constrain our parameters to satisfy 



M5 



4.29 X 10' 



(14) 



where we have used A|,(fco) = 2.43 x 10"^ for fco = 0.002 Mpc"^ [23]. In the following analysis we ensure the 
parameters satisfy this constraint. 



IV. HYBRID INFLATION: SIMULATIONS AND RESULTS 



We now report on our numerical study of the dynamics of the two fields and -0 after the end of inflation. We use 
the same rescalings as in Sec. In] With the dimensionless expansion rate H = HM~^, the equations of motion become 
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^H<l^-^^-[^^+^-^i^'<f>] (15) 

Vi + SHi/.- ^ • (16) 

The simulation space consists of a cube with comoving size L and volume V ~ discretized on a regular lattice with 
spacing Ax* = Ar (i — 1,2,3). We simulate the initial conditions of the fields as quantum perturbations around a 
homogeneous value. We take the homogeneous value of 4> to be (j)c — VX/g, the point that signals the end of inflation, 
while for ip we take ip = just as its symmetry is about to be broken. We use periodic boundary conditions. To set 
up the initial conditions, we label both free fields' normal modes by k = (27rn,j/L), where n = (n^,, rij,, rt^) and the 
Ui are integers rii = —N/2 + 1 . . . N/2. Here, N = L/Ar is the number of lattice points per side. Each free mode 
is described by a harmonic oscillator with frequency — (2sin(A;Ar/2)/Ar)^ + m^ff' where k = |k| and nies is the 
effective mass of each field given by Eq. [5] The initial conditions for the fields are then given by 



k 

^(r, i = 0) = ^ J [a.e'''-- - a^g-*''--] , (18) 

where ak is a random complex variable with phase distributed uniformly on [0, 27r) and magnitude drawn from a 
Gaussian distribution such that (|afcp) — A/2. Here we have expressed k and oj^ in units of M while Ar has units of 
M~^. The factor of A appears in our initial conditions because of our rescalings described above. The evolution of 



the expansion rate in dimensionless units is given by solving Eqs. 15 and 16 together with Friedmann's equation 



where (p) is the dimensionless spatially-averaged energy density. At each time step, we solve the three coupled 
equations for (f>, ip, and for the scale factor a(i), using H = a/a and a{t = 0) = 1. 

We discretize the equations of motion using second-order space derivatives with lattice spacing Arg — 0.05Af~^, 
second-order time derivatives with Ai = O.OIM"^, in a box with 256"^ lattice points. The effective equation of 
state parameter w = (p) / (p) is plotted in Fig. [6j Since w naturally undergoes large fluctuations, we use a window 
function to uncover its underlying behavior. As can be seen in the figure, initially the Universe undergoes a near de 
Sitter expansion with constant expansion rate ~ Hf^^ = M'*/(12Ampj) and w < —1. It then quickly transitions 
into a matter-dominated Universe which, as expected from the coherent oscillations of the inflaton about the origin, 
behaves like pressureless dust, with w = Q and H{t) cx 2/(3t). In this phase, the two fields oscillate nonlinearly 
about their respective minima with amplitudes that depend on their relative effective masses. (See, for example. 
Fig. [2] in particular the left plot) . Although we have not searched extensively, we were thus far unable to find excited 
state oscillons in a cosmological setting. The difficulty is mostly practical, since excited states require large initial 
oscillations in the </> field and this makes simulations harder to carry out. The spatially-averaged behavior of the 
two fields in a typical simulation can be seen in Fig. [7] Throughout our simulations, we set Af = 7 x 10 m-p\ and 
A = 3 X 10~^, and we investigate the evolution of the system for different values of the ratio .g^/A by changing g. At 
the same time, we adjust the value of m to be consistent with the WMAP constraint of Eq. [Til 

We identify three regimes in g^ /\. For the weakly-coupled system with g^/A ^ 1, <C and most of the energy 
of the system is transferred to 0. As a consequence, after the two fields reach the minimum of the potential, they start 
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FIG. 6: The effective equation of state parameter w = {p)/{p) for a typical simulation, showing the transition from a near de 
Sitter expansion to a matter-dominated expansion. 
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FIG. 7: Spatially-averaged evolution of both scalar fields at the end of hybrid inflation with 0c ~ M/g and t/jq ~ Mj^/X. 

oscillating with very different amplitudes. While </) oscillates with large amplitude about = 0, the oscillations of 
about its minimum at ^/^g = Mj^fX have very small amplitude. As the system evolves, loses energy to the expansion 
and is left with uniform fluctuations throughout the simulation volume. In the meantime, uniformly settles at its 
minimum and no oscillons are formed. This is the regime usually studied during tachyonic preheating [S]. 

For g^/A > 1 however, > and the amplitude of oscillations in -0 become comparable to and even larger than 
those in 0. These nonlinear, large-amplitude fluctuations lead to the formation of oscillons with an approximate radius 
of i?osc ^ hM~^ ^ a near-Gaussian profile (just as the ground state oscillons of Sec.[ll| and energy E^sc ~ 43M/A. In 
spite of the red-shifting of field modes due to the expansion, these objects remain coherent throughout our simulation 
which goes to times longer than 800M~^. This time scale should be compared with the Hubble time during their 
formation, of order ^ 200M~^. Even though tracking the oscillons for much longer times in an expanding 

Universe remains challenging numerically, our results indicate that oscillons act cosmologically as stable coherent 
objects, as has been conjectured in Ref. |13) . 
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Finally, for g'^ ^ \, cj) can be neglected in the equation of motion for and our system is equivalent to the one studied 
in Ref. ^8j, which demonstrated the existence of cosmological oscillons in a 3d, single-field model with a symmetric 
double-well potential. (Note, however, that in that previous investigation the expansion rate was kept constant for 
the duration of the simulation, while here we solve the coupled Friedmann-Klein-Gordon system self-consistently.) A 
snapshot of the energy density after oscillons have formed is shown in Fig. |8] 




FIG. 8: Energy density of the system after oscillons have formed for /X — 2. Two isosurfaces are shown, one at 10 times the 
average energy of the box (light blue) and one at 12 times the average energy (purple). The system has been evolved until the 
scale factor is a{tf) = 7. The oscillons come to rest in the comoving frame and stay localized for the duration of the simulation. 
We have used M = 1.7 x 10^** GeV, m = 1.0 x 10^^ GeV and H = 2 x 10^^ GeV. The oscillons have approximate energy 
Ease — 25 X 10^'^ GeV, confined to a radius of approximately 6 x 10"^'' cm. A visualization of a typical run can be seen at ^25]. 



Given that oscillons emerge at the end of hybrid inflation, it is important to compute their contribution to the total 
energy of the Universe. We have performed simulations until a(t) = 7 for several values of g^/A, measuring the energy 
of an oscillon by locating its core and integrating the energy around it within a sphere of radius i?shoU = 5M~^. We 
have verified that our results remain invariant using slightly higher or lower values of a(i/) and i?shoU- Since our 
initial conditions depend on random amplitudes for the fields, we average over multiple runs for each value of g^/A. 
Our results are shown in Fig. |9] These simulations have been performed using a parallel niulticore system and a 
new GPU-enhanced code, which will be described in a forthcoming publication |27j . Several features are immediately 
apparent in Fig. [oj As we remarked before, for g^/A < 1, -0 oscillates with small amplitude about its minimum, 
and thus cannot probe the nonlinear part of the potential, which precludes the formation of oscillons in these models 
[TH[T2J. Indeed, for /\ < 0.1 we do not see oscillons emerge and find that their contribution to the total energy is 
negligible. On the other hand, for g^ /\ > 3, the oscillations in (j) are too insignificant to affect the evolution of "0 and 
we find a contribution for oscillons to the total energy of about 10%, roughly independently of .g^/A. For g^/A ^ 2, 
we find a significant peak in the fraction of the total energy contributed in oscillons, with a maximum at ^ 20%. As 
we explain next, this peak is the result of parametric resonance induced by the coupling between the two fields. 



We decompose ■0 as = '>Pav(t) + S^{x,t), where ^av is the volume averaged field. Linearizing Eq. 16 with 

respect to Sip{'x.,t) and taking the Fourier transform while treating as a function of time </> = 0(t), we obtain (for 
fc > 0) 

The cosmic expansion is important only in the first few oscillations and quickly becomes insignificant during oscillon 



production. It does, however, affect the physical modes that enter Eq. 20 since they get red-shifted as fcphys = k/a{t). 
We could have taken the expansion into account by redefining x — a^ip, but since in the post-inflationary regime 
the scale factor goes as a{t) cx t^^^, the friction term due to the cosmic expansion will become negligible. As we will 
see, this choice will simplify our analysis while still incorporating the dominant effects from the expansion. 
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FIG. 9; Fraction of the energy in oscillons as a function of the couphng ratio /\. Each point represents an average over 50 
runs with error bars giving one standard deviation. 



After the fields have reached their minima, their near-sinusoidal oscillations can be approximated as ipav = "0 + 
^'cosw.^i and 4>{t) = ^cosw^t, where ■;/' ~ V'o = M/\/A and = \/2M and uj^ = g/V\M will be the fundamental 
frequencies of oscillation. Here, 4> and \E' denote the amplitude of oscillation about the minimum for each field, which 
depend on the value of /\. Inserting these into Eq. pffl we get 

5-ip{k, t) + f + a-f/i^ - 1 + eVi* cosuj^t + 3*^ cos^ uj^t + cos^ uj^t] St/jik, t) = 0, (21) 



which we rewrite as 



a{t) 



5'ilj{k,t) + (1 + hcosuj^t + ecos^ bj^t + 5cos^ to^t) 5ijj{k,t) = 0, (22) 



jI = + aV'^ -l,h = Q^i-^/ujI, e = 3^'2/w§ and 6 = For each wave number k Eq. |22 



has effectively 

only one free parameter because all quantities that appear are either constant or depend on the value of 5^ /A. 

We would like to solve the above equation analytically and look for solutions of the form i5^/'(fc,t) cx e'''=*P(t), 
where is real and positive and P{t) is a periodic function, to see the effects of parametric resonance on modes 
k/a{t). However, standard Floquet analysis I28j requires the values of 5, e and h to be much smaller than unity, 
which is not the case in this model: h, in particular, can be larger than unity. We hence solve Eq. |22| numerically, 
looking for solutions of the form Stp{k,t) cx e'^'=*P(i) and identifying the Floquet exponent /ife. We investigate values 
of < /X < 5 and 0.3M < k/a{t) < l.OM. For each value of g^/X we obtain ensemble averages for the values of 
<& and ^' from the runs of Fig. [9) which are then used to solve Eq. [22] Our results are shown in Fig. [lO] Oscillons 
of radius Rq are characterized by oscillations with wave number 2/Rq, where oscillon radii range from 2.5A/"i to 
3.5M~^ [121 [13]. And so, as previously also shown in Ref. their formation requires the excitation of modes 

with 0.57M < k/a{t) < 0.8M. In cosmological simulations, oscillons form most efficiently from initial overdensities 
of radii 4 — 6M~^, which correspond to a range of wave numbers 0.33M < k/a{t) < 0.5M. From Fig. 10 we see 
that this is precisely within the range of wave numbers that undergoes largest exponential growth, corresponding 
to the oval shaped region around g^/X ~ 2 (colored in red). These overdensities subsequently radiate their excess 
energy, shrinking slightly to condense into oscillon configurations of smaller radii, within the known range of 2.5M~^ 
to 3.5M~^ noted above. For small values of g^/X, few modes are amplified and only weakly, explaining the left part 
of Fig. [9] On the other hand, for g^/X > 3.0, modes with 0.55Af < k/a{t) < 0.8M are significantly amplified, giving 
rise to a considerable number of oscillons. Note also the near-constant bands of the amplified modes for these values 
of which clearly correlates with the flat plateau in energy in Fig. [9] (As we remarked above, for values of 

g^/X > 3.0, (j) oscillates with small and eventually negligible amplitude, while ip reaches a maximum amplitude of 
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oscillation.) Most importantly, near g^/A < 2, modes with k/a{t) ^ 0.4M undergo large amplification leading to 
copious production of oscillons. This explains the resonant peak in Fig. [9] for the contribution of oscillons to the total 
energy, -Bosc/^^tot- 
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FIG. 10: Floquet exponent /ik for solutions of Eq. 22 with different values of g /A and k/a{t) 



V. CONCLUSIONS AND OUTLOOK 



We have investigated the nonlinear dynamics of hybrid inflation cosmological models featuring two coupled scalar 
fields: the inflaton (f), with a quadratic potential, and a waterfall field ip, with a double-well potential. The fields 
couple quadratically with strength [HISlEi- Exploring the model first in static Minkowski space, we have found 
that even without cosmic expansion, the model has very rich dynamical features. We have shown that it supports 
the existence of long-lived, spherically-symmetric, coherent time-dependent structures, which we described as a new 
class of oscillons [TH [T3]. These objects appear in two possible states, which we have called ground and excited. 
In some cases, the excited state was seen to relax into the ground state, where it persisted before finally decaying. 
The lifetimes of this new class of two-field oscillons were measured to be five times longer than their single-field 
counterparts. Given that our main focus in this work is the cosmological relevance of oscillonlike structures, we have 
not provided a detailed analysis of their properties, although their enhanced lifetime is likely a consequence of resonant 
behavior. We believe they are an excellent laboratory in which to investigate the effects of resonances in coupled-field 
models able to support localized, coherent (solitonlike) time-dependent field configurations, a topic worthy of further 
study in its own right. 

Coupling the two-field model to the cosmological expansion, we have solved the Friedmann-Klcin-Gordon system 
of three coupled equations in three spatial dimensions. We started the simulation near the end of inflation, as the 
inflaton field cj) breaks the symmetry of the ip field, which relaxes to its global minimum. We note that the fields were 
initialized as Gaussian random variables, but remarkably self-organized to generate a network of long-lived oscillons 
as the fields oscillate about their minima. For certain choices of the couplings, we have shown that the long-lived 
structures persist for at least four cosmological expansion times. There were three regimes, determined by the value 
of the ratio of couplings, g^/X- The most interesting results were found for g^/X > 1, when oscillons were formed. 
Their contribution to the energy density of the Universe can be as much as 20%, in a renormalizable model with 
realistic couplings that is consistent with WMAP 7-year data. Finally, we explained the formation of oscillons and 
their substantial contribution to the total energy of the Universe as a consequence of parametric resonance. 

There are many interesting open questions to explore. Apart from studying the dynamics of the two-field model in 
static space in more detail than in Sec. ITlj it would be interesting to examine how our results would vary if we included 
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couplings to other fields in both Minkowski and expanding spacetimes. In particular, even though the oscillations of 
(j) are supposed to generate the reheating of the Universe and its transition to a radiation-dominated expansion, we 
have not incorporated the couplings to other fields that would initiate a radiation-dominated era. Hence, we found 
a transition to a matter-dominated cosmos instead. What will be the role of oscillons during reheating? If previous 
results are an indication, oscillons will act as bottlenecks for energy equipartition, delaying the approach to equilibrium 
[29] . In this case, they or other possible time-dependent coherent structures will strongly influence the value of the 
reheating temperature, providing a key constraint to our understanding of the physics of the early Universe. The 
fact that oscillons may contribute a sizable fraction of the total energy in the Universe suggests that they could have 
played a role in baryogenesis. While baryons produced before inflation are inflated away or, if produced after, are 
wiped out as the Universe equilibrates, oscillons remain out of equilibrium. It would interesting to examine if the 
time-scale associated with the oscillations of the scalar fields in oscillons, or even their lifetime itself, may be used as 
a mechanism to generate a baryon asymmetry. Work along these lines is currently in progress. 
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